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Abstract 

The magnetogenesis scenarios triggered by the early variation of the gauge couphng are 
critically analyzed. In the absence of sources, it is shown that the electric and magnetic 
power spectra can be exphcitly computed by means of electric-magnetic duality transfor- 
mations. The remnants of a pre-inflationary expansion and the reheating process break 
explicitly electric-magnetic duality by inducing Ohmic currents. The generation of large- 
scale magnetic field and the physical distinction between electric and magnetic observables 
stems, in this class of models, from the final value reached by the conductivity of the plasma 
right after inflation. Specific numerical examples are given. The physical requirements of 
viable magnetogenesis scenarios are spelled out. 
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1 Early variation of gauge couplings 



Large-scale magnetic fields were firstly detected in the optical by means of polarimetric 
observations p. The observational evidence supported the presumption that the Milky 
Way, as well other spiral galaxies, possess a magnetic field whose associated pressure turned 
out to be comparable with the cosmic ray pressure [2]. While at the beginning the opinions 
seemed to be diverging [3] the hypothesis of large-scale magnetism got progressively firmer. 
Sixty years have passed and the steady development of radio-astronomical and microwave 
observations confirm today the initial vision [U |2] which represented, incidentally, one of the 
first successful interplays between high-energy physics and astronomy. 

Since the late eighties there is mounting evidence that nearly all gravitationally bound 
systems at large scales are endowed with a corresponding magnetic field. Even if this general 
lore is not as firm as galactic magnetism, the general evidence of magnetic fields associated 
with Abell clusters is becoming progressively more compelling (see [1] for a general review). 
The existence of large-scale magnetic fields in super-clusters is, at the moment, not yet 
established because of intrinsic observational drawbacks which might become less severe with 
the advent of a new generation of experiments such as the ambitious Square Kilometer Array 
(SKA) [5]. If confirmed, the recent Auger observations [7] might constrain the magnetic fields 
present in a cocoon of 60-70 Mpc. 

One of the intriguing questions present from the very beginning [2] concerns the origin 
of large-scale magnetism. This question became more acute since all spirals have magnetic 
fields with comparable strengths and similar correlation scales [8]. There are certainly im- 
portant observational issues related to the latter problem and at least some of them could be 
addressed in the framework of the SKA project (see, for instance, third paper in [S]). At the 
same time only galactic observations will not tell us, with reasonable confidence, if large-scale 
magnetism is a phenomenon which is primordial or recent. In the present context primordial 
means that the magnetic field had to be already present before photon decoupling. The true 
problem before us today is, therefore, not the endorsement of a particular model. The true 
issue before us today concerns the presence rather than the absence of large-scale magnetic 
fields before photon decoupling [6j. Are we able to exclude, to a certain confidence level, 
that pre-decoupling magnetic fields did distort the observed Cosmic Microwave Background 
(CMB) anisotropics and polarization? This question motivated a bottom- up approach to 
magnetized CMB observables. Recent results show that large-scale magnetic fields larger 
than 3.5 nG are excluded to 95 % C.L. and for magnetic spectral indices = ^-^-1,1 (see 
[6j and references therein). While further studies along this direction might either confirm 
or infirm the obtained results, it is important to stress that, absent specific bounds and 
observations on pre-decoupling magnetism, the primordial nature of large-scale magnetic 
field will just become the arena of theoretical discussions. Having said this, it is interesting 
to speculate on the origin of the pre-decoupling magnetism even if the final referee, in this 
dispute, will not be theoretical opinions but the observational evidence. 
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There are, indeed, various models for the origin of large-scale magnetism based on diverse 
physical pictures j^. The mechanisms based on the parametric amplification of minute 
fluctuations of a putative Abelian field allow to obtain correlation scales which are naturally 
large. In the present paper the attention shall be confined on those magnetogenesis scenarios 
which are based, in a way or in another, on the assumption that the hypermagnetic and 
hyperelectric fields are coupled to a spectator field whose evolution is dictated by a quasi-de 
Sitter dynamics [9j. Recently various aspects of these models have been investigated [10]. A 
related class of models stipulates that the the gauge coupling takes part in the dynamics of 
the background and this can arise either in the context of conventional inflationary models 
[TT] or in the context of string cosmological scenarios [T2| . Only some of the considerations 
developed here are apphcable to both classes of models and probably only the results of 
section [2] are general enough to embrace all the different possibilities. 

Consider, for simplicity, a spatially fiat Friedmann- Robert son- Walker geometry whose 
line element can be written, in the conformal time parametrization as 

ds^ = G^^dx^dx'' = a\T)[dr^ - dx% G^^ = a\T)r]^^ (1.1) 

where r denotes the conformal time coordinate and r]^,^ is the Minkowsky metric. In what 
follows Vq will denote the covariant derivative with respect to the metric (11. ip . In the 
absence of any concentration of particles charged under a putative (Abelian) gauge field the 
action can be written as 
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Sy = -— / d^xV^ 
lovr 



A(^)r,/3F"^ + A(V')r„^r"^ 



:i.2) 



where both, X{ip) and X^tp) depend explicitly on the field ip whose action and potential can 
be specified in a given scenario; furthermore, in Eq. (II. 2p 



a/3po- 

Y^p = V[„y;3] = d^^Y^], = i^^Y,,. (1.3) 

The explicit form of the field strengths in terms of the electric and magnetic fields are given, 
respectively, by Y^i = a? Cj and by l^j = —a? e-ijk while the dynamical equations and the 
Bianchi identities are 

v« (^Ay"^ + Ay"^) = 0, v„y"^ = o. (1.4) 

In terms of the canonical hyperelectric and hypermagnetic fields the action becomes: 

Sy = I f d^x[E^ - B% B = J^a^h, E = J ^ a" e; (1.5) 
2 J V 47r V 47r 

the coupling X{ip) will now be set to zero even if, as we shall remark, the considerations 
on the duality of the electromagnetic spectra obtained from the parametric amplification of 
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vacuum fluctuations can be generalized to the case where A 7^ 0. In the absence of charged 
particles the evolution equations of E and B can be written as 

Vx(Va5)-A(VaJ)^0. |(^)Wx(|)^0. (1.6) 

Under electric-magnetic duality, i.e. 

E^-B, B^E, (1-7) 

the first of Eq. (11.61) goes into the second (and, vice- versa, the second equation is transformed 
into the first one). In this sense plays here the role of inverse coupling 5fy/(47r) = 
tty = 1/VX. Note that the action (11.61) is not invariant under (II. 7p (see also [13] for 
detailed discussions on the canonical formalism applied to duality). In the present context 
the transformation of Eq. (I1.7p will simply be used to compute the electric and magnetic 
power spectra when the gauge coupling evolves with the conformal time. If A 7^ in Eq. 
(II. 2p the resulting theory can be shown to possess a continuous SL{2, R) symmetry which 
generalizes ( 11. 7p . Some of the considerations illustrated here can be generalized to the full 
SL{2, R) case but this will not be our primary purpose. 

The plan of the present paper is the following. In section [2] it will be shown how the 
electric and magnetic power spectra can be related by duality transformations both exactly 
and approximately (i.e. when the relevant wavelengths are larger than the Hubble radius). 
The analysis will be conducted primarily in the case of single field inflationary models even 
if occasional remarks on different models will be made. In section [3] it will be shown how the 
finite value of the conductivity breaks the duality symmetry and distinguishes the magnetic 
degrees of freedom from the electric ones. In a more formal sense the presence of the 
conductivity induces an effective decoherence of the amplified electromagnetic fluctuations. 
In section H] numerical examples will be discussed and the parameter space of the model will 
be be studied. The concluding remarks are collected in section [5l 



2 Electric-magnetic duality and power spectra 

In time-dependent (and conformally flat) backgrounds such as the one of Eq. (II. ip the 
quantization can be easily performed, for the present ends, in the Coulomb gauge. The 
Coulomb gauge condition (i.e. 10 = and V ■ F = 0) is preserved under a conformal 
rescaling of the metric (11.10 : the Lorentz gauge condition does not have the same property 
[H]. Furthermore, in the Coulomb gauge, the duality properties are easily discussed □ [13]. 

^It is of course possible to keep also Yq and the longitudinal part of the Abelian vector potential, i.e. 
Il- It can be easily shown, however, that the transverse variables are decoupled from the longitudinal and 
from the gauge contributions. Indeed, without any specific gauge fixing, the extremization of the action with 
respect to Yq implies that = VIq. 
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The action of Eq. fll.2p can then be written as 

Sy = J dTLyir), Ly{t) = j d^x Cy{x,t) 



(2.1) 



where, in the case A = 0, 



(2.2) 



and y = yX/{ATT)Y. The canonical momentum conjugate to y can be easily obtained from 
Eq. (12. 2p and it coincides, up to a sign, with the canonical electric field, i.e. 



n = y -^y = -E 



(2.3) 



while the relation of y to B is simply given hj B = V x y. The canonical Hamiltonian is 
then given by 



Hy(t) = - I d^x 



7f^ + 2 



The Fourier mode expansion for the canonical fields 



TT-y + diy-d'y 



(27r)3/2 

can be inserted into Eq. (12. 4p and the resulting expression is 



d^k vAt) e 



-ik-x 



(2.4) 



(2.5) 



Hy{t) = ^ I d'k 



■ y^k + 7r„fe ■ yg) + k% ■ y^ 



(2.6) 



Equation (12.61) is invariant under the transformation 

1 ^ 



fc' yk j^^fc' 



1 

k' 



-ky_ 



fe' y~k ^ U^-k^ 



(2.7) 



where k = \k\. The Hamilton equations derived from Eq. (12.61) become: 



-yh 



(Va)' 



(2.8) 
(2.9) 



By applying the transformation (12. 7p . Eq. (12. 8p is transformed into Eq. (12.90 and vice versa. 

In cosmological backgrounds, any canonical transformation depending explicitly upon 
time does modify the form of the Hamiltonian without altering the Hamilton equations. 
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While such an occurrence is harmless at the level of the evolution equations, the different 
forms of the Hamiltonian can make the difference when setting initial conditions at a finite 
time during the inflationary epoch. A similar problem arises for spin-0 [15] and spin-2 
fluctuations of the geometry [16] and it is related to the so-called transplankian ambiguity of 
inflationary initial conditions. In the present context pre-inflationary initial conditions can 
induce explicit sources which break the symmetry of the system of equations given in Eqs. 
^Ml and ([23]). 

The canonical flelds can be promoted to quantum operators (i.e. yi jji and tTj — > vrj) 
obeying (equal time) commutation relation^: 



[yi{Xi,T),7Cj{x2,T)] =iAij{xi -X2), 



AiJxi - X2) 



(2^ 



e'''<"^^P,^{k), (2.10) 



where Pij{k) = {6ij — kikj/k"^). The function Ajj(xi — X2) is the transverse generalization 
of the Dirac delta function and such an extension is mandatory since V ■ E = because of 
the Gauss constraint and V ■ y = because of the gauge conditioE0. The fleld operators will 
then become 



7li[X,T) 



d^k 



(27r)3/2 ^ 



t ik-x 



d^k 



i:^t\k) 



(27r)3/2 

where, from Eq. (12.91) . f\ and Qk obey: 

Jk — 9k-\ y^'''" 



fkir)dk,ae *'-" + /:(r)aI„e 



^ f, ^,—ik-x I ^ik-x 

gk{r)ak,ae +9kV^)0'ha(i 



(2.11) 
(2.12) 



9k 



-k'fk - 



9k- 



(2.13) 



The mode functions /a:(t) and gk{T) must also satisfy (Wronskian) normalization condition 
which follows from enforcing the canonical commutators: 



fk{r)9l{r) - f;iT)g^ 



ikir) = I. 



(2.14) 



In Eqs. f l2JT]) and (I2J2D ef\k) must obey Ea (iT\k) e'f'ik) = Pij{k). Finally the creation 
and annihilation operators appearing in Eqs. (12.111) and (12.121) obey [a^„, dt^] = 6ai3S^^\k + 
p). The state annihilated by dk^a, i-e. |0), is indeed the vacuum but only in the limit 
T —00. Such a limit is, strictly speaking, not physical since expanding de Sitter space 
is not geodesically complete. This consideration implies that the normalization procedure 
must be conducted more carefully whenever the number of inflationary e-folds is close to 
minimal. Corrections depending upon the flnite normalization time can appear in the power 
spectrum [T5l [T6] . These corrections will be neglected here but can be important in more 



Units h = c = Kb = ^ will be assumed hereunder. 
"^Notice that A and A are here assumed to depend only upon the conformal time coordinate r. 
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refined treatments. The two point function of tlie electric and magnetic fields for equal times 
but different spatial locations is [9j: 

(0|5,(f,r)S,(f + f,r)|0) = j dlnkVB{k,T)P,^{k)^—, (2.15) 

^ ^ p Sill hjT" 

{0\Eiix,T)Ej{x + r,T)\0) = J dlnkVEik,T)P,j{k)^—, (2.16) 
where VB{k,T) and VE{k,T) are, by definition, the magnetic and electric power spectra: 

VE{k,r) = |^|/.(r)P, Vs{k,T) = ^,\9k{r)\'. (2.17) 

The expectation values of the energy density as well as of all the other components of the 
energy-momentum tensor can be derived in similar ways. The energy momentum tensor can 
be written directly in terms of the canonical fields E and B introduced in Eq. (11.51) : 

r„° = ^(B'' + B'), T-^^iExB)' (2.18) 

By averaging the different components of the energy- momentum tensor we get: and it is 
given by: 



(0|T°|0) = py(r) = l|rflnfc 



rE{k,T)+VB{k,T) 



{0\Ti\0) = -py{T)6i, py = ^, (2.20) 



where it has been used that, for a generic function h{k, r) 

4:71 



J (fk kik^h{k, r) = -^6i J k^ h{k, t) dink, J d^kk'h{k, t) = 0. (2.21) 



When the gauge coupling evolves in time as 



ay(r) = «y(r,)(--j , gy(r) = gyiT,)^--J , (2.22) 



the solution of Eq. (12.131) can be written as 



Mr) = V^Hi'\-kr), N = y|e-(^+V2)/2^ (2.23) 



9k{r) = -Ar^V^Hl'\{-kT). (2.24) 

Equations (I2.23P and (I2.24p obey the Wronskian normalization condition of Eq. (12.141) as 
it can be easily verified by recalling the properties of the Hankel functions [20l [21] . In the 
hmit r ^ -oo Eqs. (^2E^-^M) become 

hm /,(r) = ^e-^^ \\m g,{T) = -tJ '^e-''\ (2.25) 
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If u = 1/2, A is constant and the asymptotic solutions of Eq. (I2.25P are actually exact 
solutions. When z/ = /^(t) diverges logarithmically in the limit of small argument, i.e. 
|fcr| -C 1. Using Eqs. fl2.23p and fl2.24p the electric and magnetic power spectra of Eq. (12.170 
become 

V^{k,r) = f{-kT)\H;^},{-kT)\', r^{k,r) = ^{-kr)\HW(-kr)\'. (2.26) 
Svr Sit 

Let us now perform a duality transformation and show that the electric and the magnetic 
power spectra are interchanged while the energy density is invariant. Consider, therefore, 
the following transformation 

grir) ^ -grir) = ~gy{r) = ~gyin)( . (2.27) 

9Y{r) \ T-J 

By repeating the same steps discussed above with the transformed gauge coupling it is clear 
that 

Mk,r) = f{-kr)\Hl^\{-kT)\\ ^k^r) = f{-kr)\Hi'\~kT)\' . (2.28) 

But Eqs. (12:221) and (12:28|) imply that 

i) = l-u, VE{k,T) =VB{k,T), VB{k,T) =VE{k,T). (2.29) 

The result of Eqs. (12. 28^ and (12.29^ follow immediately by appreciating that [20l [21] 

\H!^\i-kr)\' = \H['li-kr)\\ H^^liz) = e^^^ H^^\z) (2.30) 

where /i and z denote, respectively, a generic (real) index and a generic argument of the 
Hankel function. The second relation of Eq. (I2.30p simply states that a generic Hankel 
function of first kind with index — /x equals the corresponding Hankel function of index +/i 
up to a yU-dependent phase which is immaterial as far as the square moduli are concerned. 
For the Hankel function of second kind the same kind of relation holds and it follows by 
taking the complex conjugate of the second relation in Eq. (I2.30p . On the basis of the 
obtained results it is also clear that under Eq. (I2.27P 



Py{t) pY{r) = ^ / c^lnA; 



VEik,T)+VBik,T) 



pY{r) 



Py PY{r) = Py{t). (2.31) 

It is useful to define, for immediate peruse, the spectral energy density in units of the critical 
energy density 



PE{k,T)+PB{k,T) 



(2.32) 
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where pcrit = 3H'^/{8nG) is the critical energy density. Using the resuhs of Eq. fl2.3ip . 

Consider now the case when the evolution of the geometry evolves according to a quasi-de 
Sitter stage of expansion driven by a (single) scalar degree of freedom tp. After the quasi-de 
Sitter stage the conventional lore suggests that the radiation dominated stage of expansion 
takes place almost suddenly. In what follows the latter possibility will be examined together 
with different physical scenarios like the ones stipulating that the transition to radiation has 
a finite duration (see also section [3]). It will be also worth examining, for the present ends, 
all the cases when the transition to radiation is delayed by a long phase where the speed of 
sound of the plasma is either softer of stiffer than radiation (see also section Hj). 

During the quasi-de Sitter stage of expansion the relevant slow-roll parameters are given 

by 

^ = -i « ^- " = w ^< 

Denoting with V{ip) the infiaton potential, the slow-roll equations can be written in the form 

H^Ml ~ —V, SHip + ~ 0. (2.34) 
3 dip 

Using Eqs. (12.331) and (12.341) it is rather standard to derive the following relations 



yj. V = e-% V = M;^. (2.35) 

where the reduced Planck mass Mp = Mp / v^Stt has been introduced. To make the forth- 
coming statements quantitative (rather than, as sometimes done, just qualitative) we recall 
that the power spectra of scalar and tensor fluctuations are assigned as 

Vn{k) = A^^(^-j , Vt = rTAn[j^) , (2.36) 

where An and r^Aii are, respectively, the amplitudes of the scalar and of the tensor modes 
at the pivot scale kp = 0.002 Mpc~^ and where the spectral indices are related to e and fj as 
well as to rx 

nT = -2e, n, = l-6e + 2r], n = 16e = -Sut. (2.37) 



The WMAP Syr data alone in the light of the vanilla ACDM scenario imply [T8j A-r = 
(2.41 ± 0.11) X 10~^ and provide an upper limit on tt and e, i.e. tt < 0.58 and e < 0.036 
(95 % CL). The bounds on e can be improved by combining the WMAP5yr data with other 
data sets. For the present purposes, explicit numerical estimates will assume, as it will be 
seen. An = 2.41 x 10"^ and e = 0.02 implying 

l-SlxlO-^^fTV^V ,,,-^^,,., ). (2.38) 



■ V0.02/ V2.41 X 10-9 
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The slow-roll parameters enter the relation between the conformal time coordinate r and 
the Hubble rate. More precisely, recall that 



dt 
a(t) 



aH 



+ e 



da 



(2.39) 



where the second equality follows after integration by parts assuming that e is constant (as 
it happens in the case when the potential, at least locally, can be approximated with a 
monomial in (f). Since / dt/a = J da/{a'^H), Eq. fl2.39p allows to express aH in terms of r 
and e, i.e. aH = — l/[r(l — e)]. Using the latter relation Eq. (12.321) can be expressed as 



H' 



X 



3M^ (1 -e) 



X 



1 - e 



H 



(1) 



X 



1 



where e is the slow-roll parameter introduced in Eq. (I2.33p and the notation 

k X 



kr 



aHil 



1 



(2.40) 



(2.41) 



will now be employed. An interesting limiting case is the one where x/{l — e) -C 1 where 
the explicit expression of f2y(x, u) can be evaluated by using the small argument limit of the 
Hankel functions. This limit corresponds to physical wavelengths which are larger than the 
Hubble radius during the quasi-de Sitter stage of expansion and in the absence of Ohmic 
currents. For practical purposes it is useful to define the function 



C{z,e) 



37r2 



'1 _ e)2^-i22^ 



(2.42) 



where z is a real number which is otherwise unrestricted. If > 1/2 the spectral energy 
density will be given, respectively, by 



Q'^\x,u) 



X, u 



El 

Ml 



Ciu, e)x 



5-2u 



Civ, e) X 



b~2u 



1 + 
1 + 



r^(z/-l)a:^ - 
4r2(z/)(l -e)2. 
r^fl - v) 22-4- 



z/ > 1, 



Au-1 



X 



iv-2 



(2.43) 



- < z/ < 1.(2.44) 



The case of increasing gauge coupling can be immediately obtained from Eqs. ( 12.43^ and 
(I2.44p . Indeed, recalling the transformation of Eq. (12.271) and bearing in mind the relation 
z> = 1 — z/ we shall have that 



H^ 



3+2v 



1 + 



-Z/ X 



4r2(l-z/)(l-e) 



z/ < 0, 



(2.45) 



^C(l-z.,e)x- 



1 + 



r2fz/)2 



X 



2-4u 



THl-u) (1-e 



2-4u 



0<z/<-,(2.46) 



where the cases z/ < and < z^ < 1/2 correspond, respectively, to z/ > 1 and to 1/2 < z/ < 1. 
It is useful to notice that the duality of the spectra can be productively used to relate the 
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four cases discussed above which reduce effectively to two cases since, as it can be easily 
seen, 



n'^^ (x, u) = n'y' (x, 1 - z/) , n'y' {x, u) = n^y> (x, i - u) 



id), 



The transformation z/ — > z/ = 1 — z/ implies that in the cases z/ 
energy density is the same and it is given by: 



and u 



flyix, Z/) 



El 

Ml 



a;=^C(l, 



1 + 



X 



(1 



In^ 



X 



(2.47) 
the spectral 

(2.48) 



Concerning Eqs. f l2.43l) -( !2.44l) and Eqs. f l2.45l) -( l2.46p few comments are in order. It can be 
easily verified that the spectral energy density is increasing for —3/2 < z/ < 5/2 while it is 
decreasing for v < —3/2 and for v > 5/2. The case of flat spectrum can be realized both for 
z/ = — 3/2 as well as for z/ = 5/2; finally the closure bound seems to imply, grossly speaking, 
that g g 

iy>---uj, v<- + uj, cj~0.3. (2.49) 

The estimate of Eq. fl2.49p follows by recalling that when the spectra are increasing the 
most relevant constraint stems from Xmax — 1; when the spectra are decreasing the critical 
density bound should be applied for x = Xmin where 



3.149 X 10 



-30 / 



(2.50) 



The estimate of Eq. (12.501) takes into account an intermediate (expanding) phase character- 
ized hj n = 2/[3{wt + 1)] where wt is the barotropic index characterizing the intermediate 
phase while H,. marks the onset of the radiation-dominated stage of expansion. In the limit of 
sudden reheating Hj. = H. By neglecting all the numerical pre-factors in Eqs. fl2.43l) - fl2.44l) 
and Eqs. (I2.45p - (l2.46p the value of u) can be estimated as: 

log (ff/Mp) 



UJ 



log x„ 



(2.51) 



By taking into account all the pre-factors the estimate uo ~ 0.23 for e = 0.02. Note that the 
estimate of Eq. (12.511) is based, just for illustration, on the sudden reheating approximation, 
i.e. = H in Eq. (12.501) . The limits stemming from Eqs. (12.511) (as well as from Eq. 
f l2.49p ) just mean that, as soon as u gets larger than —(3/2 + u) (or, equivalently, smaller 
than 5/2 + u) the dynamics of the gauge coupling should be discussed self-consistently by 
using the following set of equations 



Mp{n'^-n') 
if" + 2n^' 



2^^ 



1 

^3 



1/2 



+ PY, 



ij" + 271^' + 



dV_ 

dip 
dW 

dil> 



0, 



a? + 



1 d\n\ 

2a2 dil) 



(2.52) 
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where, as in Eq. fl2.35p . the reduced Planck mass Mp = I/vSttG = Mp/yoTT has been 
introduced. In Eq. (12.521) W{ip) denotes the potential of ip, V{ip) (as in Eq. (12.350 ) is the 
inflaton potential while Ti = a'/ a. 

The electric and magnetic spectra obtained in the present section assume that the system 
of governing equations are invariant under the electromagnetic duality. However, in the 
transition regime from inflation to radiation charged particles are created. Therefore the 
conductivity switches on from zero to a finite (large) value and the duality symmetry is 
explicitly broken. 

The analysis leading to Eqs. (I2.43p -( l2l44ll and to Eqs. (I2.45l) -( l2^ is purely kinemati- 
cal.This means that a putative form of the gauge coupling has been assumed by completely 
disregarding how the dynamical scenario could be realized. Any limit derived on the spec- 
tral index just highhghts the necessity of a more complete description which was attempted, 
incidentally, in [9] (see also [T7j for the analysis of electromagnetic backreaction effects in a 
different but related context). Conversely, the (qualitative) conclusions of [22] only demon- 
strate the inadequacy of an approximation for a limited range of the parameter space of the 
toy model under discussion. 

In the case of pre-big bang models [12] the situation is different from what has been 
described so far and the bounds discussed above do not arise. For instance, in the toy 
model used in [17] the evolution of the gauge coupling and of the geometry is regular. The 
asymptotic evolution can be written, in the Einstein frame, as: 




a(r) ~ a_ J , a_ = e 



2(73-1-1) 



ln2-V31n(^)-y31n(-^) 



= f , = - — . (2.53) 

Zt t 

for r — s> — oo. The gauge coupling will then be given, in this context, by e*^"/^ which implies 

u = (1 — \/3)/2 in the equations for fiy(r, u). In this case the critical density bound is never 

violated since, to leading order in x < 1 

ny{x, u) ^ ^x'-^\ ^ = 4^2 « 1- (2-54) 
Mp Mp a^Mp 

In this case the result seems to be different from the qualitative estimates of [22j and in 
agreement with former results [12]. Similar estimates can be carried on in the so-called 
string phase of the model [12]. Back-reaction effects of the gauge fields can be relevant, in 
bouncing models of pre-big bang type, but in the opposite regime. The amplified quantum 
fluctuations can re-enter the Hubble radius when the geometry is still dominated by the 
dilaton. In this case the energy density of the amplified fiuctuations decreases more slowly 
than the background and the result can be a kind of gravitational reheating which has been 
studied, in detail, in [T7] . 
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3 Breaking of the duality symmetry 



The transition from inflation to radiation can be smoothly parametrized in terms of the scale 
factor of a conformally flat Friedmann-Robertson- Walker background 



a(r) = r + + rf. (3.1) 

From Eq. (13. ip it follows that 

a(r):.2(-^), Ito a(r) (1), (3.2) 

implying that for r <^ — ti the geometry is of quasi-de Sitter type while for r ^ ri the 
geometry is dominated by radiation. The transition regime is associated with the scale 
Ti. The smoothness of the transition is necessary if the conductivity has to be included 
consistently. In the presence of an Ohmic current the evolution equations discussed in Eqs. 
(I1.6P are modified and they become @ 

V X {^/XB) = ^[v^^] + 4vrJ, (3.3) 



d_ 



B 

71 



V X 



E ■ 

71 



0, (3.4) 



where J is the Ohmic current arising because of the interaction of the electromagnetic 
field with the ambient plasma. As pointed out in |2| the hypercharge field projects on 
the electromagnetic field as Af^ = cos6'w3^i- In what follows, however, we will simply 
assume the presence of a relativistic electromagnetic plasma with positively and negatively 
charged species whose masses are much smaller than the corresponding kinetic temperatures. 
Defining as n± = a^n± the comoving concentrations of positively and negatively charged 
species we shall assume that the Universe during and after inflation is globally neutral, i.e. 
n+ = n_ implying that V ■ E = 0. The Ohmic current can be written, according to the 
conventions adopted above, 

J = a^j = a^acS = —/=E, a = a^a. (3.5) 
V A 

Note that a = acxc is the rescaled conductivity which can be computed for a relativistic 
plasma with evolving gauge coupling. It is also relevant to remark that J scales with 1 / V^: 
this occurrence can be intuitively understood by recalling that, when the gauge coupling is 
not dynamical, the current is proportional to the electric charge, i.e. 1/VX in our notations. 

The conductivity is vanishingly small during inflation and it jumps to a finite value after 
the inflation. Similarly the gauge coupling will evolve dynamically during inflation and it 



^For practical reasons we shall use hereby Gaussian units with the An factor in front of the current. 
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will eventually freeze during the subsequent phases. An interpolating expression for A(r) 
can be written as 

where a = (1 — 2i^)/3. From the algebraic form of Eq. (13.61) it follows that 

hm A(^) = ^ 



hm A(^) = - hm A(^) = (3.7) 



In the same interval of (rescaled) conformal time the conductivity must jump from zero to 
a finite value which coincides with the temperature and this behaviour can be modeled as 

ae(.) = ^^(z), e{z) = U^ + ^^X. (3.8) 



In Eq. (13.81) T^h is the reheating temperature. By increasing (3 the transition becomes 
sharper. Indeed as (3 increases (Tc vanishes more accurately during the inflationary stage of 
expansion. 

In the presence of the conductivity the evolution equations of the mode functions are 
modifled as follows 

, ,2. 47r(T (v^)' 
9k = -k fk yQk '^3^' 

fk=^-^fk + 9k. (3.10) 

The evolution equation of the mode functions in the conformal time coordinate r becomes 
simply: 



( ■/X)" ^ im ( -/X)' 



When we pass from t to z = t/ti the terms containing cr, for power counting reasons get 
multiplied by ti = l/7ii = ^/{cliHi). The relevant combination can then be written 

inari _ A7rac{z)Tia{z) _ Wtc'^Ta a{z)e{z) 

A ~ Xiz) ~ Hi e^X{z) • ^ ' 

where we used the fact that a{z) = ac{z)a{z) and that 

, 7]-2 / 90 \ ^Z"' / — 

HfMp = —N,sTX, = — — ^/HiMp, (3.13) 



90 VTT'^A'efr 

where N^s is the effective number of relativistic species which will be assumed to coincide, 
for simplicity, with 106.75 as in the standard electroweak model. 

Conductivity effects might also play a role at the onset of the inflationary evolution. The 
maximal number of inflationary e-folds accessible to our observation correspond to the min- 
imal number of e-folds required to relax the problems of the standard cosmological model. 
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If the case of a standard post-inflationary thermal history A^max = 67.951 + 0.25 In (ttc^t^). 
For our fiducial set of parameters iVmax = 62.3. This estimate is consistent with the result 
obtained in [19] for the number of e-foldings before the end of inflation at which observable 
perturbations were generated. A deviation from the standard post-inflationary history (with 
the plausible addition of a stiff epoch) can reduce iVmax even by 10 e-folds. If the number 
of inflationary e-folds is (9(63) the initial state of the electromagnetic perturbations is sen- 
sitive to the pre-inflationary expansion which is unknown and which could have contained a 
relativistic plasma. In the latter case the electric degrees of freedom are already suppressed 
when the initial conditions of the mode functions are set. 

4 Phenomenological considerations 

In the case of conventional inflationary models the evolution of the geometry starts, typically, 
from a highly curved stage of expansion. If the duration of inflation is not minimal the Hubble 
rate at the onset of inflation can be close to the Planck scale. Assuming e = 0.02 from the 
observational data and supposing, for instance, that inflation lasted a number of e-folds six 
times larger than the minimal amount (i.e., overall, A^efoids ~ 360 or even larger) the Hubble 
rate at the onset of inflation could have been O(10~^'^)Mp or larger. Since the inflationary 
evolution commences in a regime of strong gravitational coupling, it is not unreasonable that 
also the gauge coupling could be strong at the onset of the dynamical evolution. Consider, 
as an example, the case a = —4/3 and u = 5/2. The evolution equations of the electric 
fields can be integrated numerically by imposing quantum mechanical initial conditions at 
the time ri. The choice of quantum mechanical initial conditions is appropriate due to the 
largeness of A^'efoids- On the contrary, if the pre-inflationary evolution is dominated by a 
relativistic plasma the minimal number of e-folds would not guarantee the plausibility of 
quantum mechanical initial conditions at T[. It is now practical to introduce the rescaled 
wavenumber k = kri where Ti sets the scale of the transition from the inflationary phase to 
the radiation dominated phase (see Eq. (13.11) ). At the onset of the numerical integration all 
the wavelengths are all inside the Hubble radius, i.e. nzi > 1 where Zi = t^/ti. In Fig. [T]the 
result of the numerical integration is illustrated for the case z/ = 5/2. In the plot at the left 
the conductivity vanishes and, as expected, the spectrum of the electric field scales according 
to the analytical result, i.e. V-e — i^"^- In the plot at the right the conductivity is included 
and the electric power spectrum is exponentially suppressed. The exponential suppression 
does depend upon the evolution of the conductivity across the matter radiation transition. 
In Fig. [1] the plot at the right assumes the parametrization of Eq. (13. 8p with /? ^ 1 (in 
practice (3 = 8). The latter choice guarantees that during the inflationary phase a vanishes 
rather accurately. In the opposite case (i.e. P ~ the presence of the conductivity can 

modify both the amplitude and the slope of the spectra. In the simplified model discussed 
here the critical requirement is therefore that a/X must vanish much faster than, say. A"/ A 
for T <^ — Ti. In Fig. [2l Eq. fl2.14p is numerically illustrated. The result of Fig. [2] implies 
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Figure 1: The electric power spectrum with (i.e. o" ^ 0) and without (i.e. a = 0) conductiv- 
ity. 

that the evolution of the canonical operators is such that [iii.T^j] — > for r > ri. Because 
of the presence of the conductivity the Wronskian is driven to zero. This means that out of 
the two solutions of the system only one survives, i.e. the one related to the magnetic part. 
The vanishing of the Wronskian signals the transition to the classical dynamics where the 
magnetic field operators become Gaussian random fields, i.e., 

B^{x,t) = j d'kBAr), {BAT)B,{p,r)) = '^VB{k) P.,{k) 6^'\k + p), 

(4.1) 

In Fig. [3] the spectra of the magnetic fields are illustrated always in the case of decreasing 



V = 5/2, A„ = 2.41 x^Q~^,altO 

H 




t/t^ 

Figure 2: The vanishing of the Wronskian is illustrated. 
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gauge coupling. It is practical to assign the magnetic power spectra with the same con- 
ventions employed in the case of the power spectrum of curvature perturbations (see Eq. 
dlSSD): 

k 



n.R-l 



(4.2) 



where /cl = 1 Mpc^^ is the magnetic pivot scale and is the spectral amplitude at ki^. 
In Fig. [3] the cases tt-b = 1 and n-^ = 1.5 are explicitly illustrated. When the transition to 
radiation is delayed (in the case of a small (3) the magnetic fields are more suppressed than 
in the case of (3 ^ 1 (which is the one illustrated both in Figs. [T]and[3])). The results are 
most easily normalized in terms of where Hi = 7^(0) since Ti-ir) = 1 / -y/r^ + rf . In the 
case of flat spectrum the amplitude ^b is solely controlled by the amplitude of the curvature 
perturbations and it can be written as 



A, 



P7(^o) 



18.49 e^Te 



2.725 K 



p^(ro) = 2.29 /iG 



(4.3) 



where ^^(to) = is the energy density of the CMB photons. Equation (14.31) implies 

that v^v4b = 6.840 x 10~^^ G = 0.068 nG which agrees with former estimates [9]. To convert 
the various results in physical units it is useful to recall, as mentioned in Eq. (14. 3p . that the 
energy density of the CMB photons does correspond to a putative magnetic field of about 
3 fiG (1 fiG = 10~^ G) which is, in turn, comparable with the energy density of the cosmic 
rays as already argued many years ago[2j. In the case nB > 1 the analog of Eq. (14. 3 p can be 



n„=1, e = 0.02, A„ = 2.41x 10""^, a^Q 



n„= 1 .5, e= 0.02, = 2.41 x 1 0"", a*0 




O 




K= 10"' 



■ ■ K= 10"' 



- K=10"' 



X/X, 
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t/t. 
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Figure 3: The magnetic field spectra for two different values of the spectral index (i.e. wb = 1 
and riB = 1-5) and for different wavenumbers. 

written as 

^/ \(5-"b)/4 X (l-2M)(nB-5)/2 
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where Hi denotes the Hubble rate at the end of inflation, Hj. is the Hubble rate at the onset 
of the radiation-dominated phase and n parametrizes the (average) exponent of the scale 
factor between Hi and H,-. The function F{nB, e) is given by: 

27.362 - ub [26.307 + 0.434 ln(l - e)] + 2.171 In (1 - e) + 0.868 In [1(3 - nB/2)]. (4.5) 

It can be readily checked that in the case Hi = H,- Eqs. (14.41) and (14.51) reproduce the result 
of Eq. (14.31) when ne = 1 and for the typical values of An and e. Equations (14. 4p and (14.51) 
hold, both in the case of decreasing and in the case of increasing gauge coupling. However 
the relation of u to ris changes slightly in the different cases. If > 1/2 we shall have 
u = (3 — ^6/2) with ne < 5. If < 1/2 we have to distinguish two cases. If0<z^<l/2we 
shall have that u = {3 — 71^/2) (with 5 < tib < 6); if z/ < then v = {ub — 6)/2. 

If the duration of inflation is close to minimal the initial conditions of the electromagnetic 
fluctuations depend upon the pre-inflationary evolution much more than the initial conditions 
of, for instance, the scalar and the tensor modes of the geometry. Suppose, for instance, 
that the pre-inflationary expansion rate was dominated by a relativistic plasma. Suppose, 
furthermore, that inflation lasts around 63 e-folds implying that, at the onset of inflation, the 
Hubble rate was indeed O{10^^). By assuming a sudden transition from the pre-inflationary 
regime to the inflationary phase, we can estimate the upper limit on the temperature which 
is given by ~ QOif^Mp/ (Tr^A^^eff) where A^cfr is the effective number of relativistic species. 
Since the gauge coupling is small, no matter how small the temperature is, the conductivity 
will be very large since, approximately, as we saw ac — T/ay where ay = gy/i^Ti). In this 
situation the electric fields will be suppressed in comparison with the magnetic fields and 
the initial conditions for the mode functions will have to reflect this occurrence. 

It is now the moment of assessing the orders of magnitude of the magnetic fields obtained 
within this class of mechanisms. Let us first consider the case of an exactly scale-invariant 
spectrum (i.e. ub = 1) and let us measure ^^Ab in units of nG. The amplitude of the 
scale-invariant spectrum does depend, primarily, upon e and upon An- The two horizontal 
lines both in the left and in the right plot of Fig. H] illustrate the minimal case where 
( = Hj-fHi = 1. In this situation the Hubble rate at the onset of the radiation epoch 
coincides with the Hubble rate at the end of inflation. It can happen that the onset of the 
radiation-dominated epoch is delayed either because the reheating process is prolonged or 
because of the presence of an intermediate phase. This situation is parametrized, in Fig. H] 
by the value of C < 1 and by the barotropic index during the intermediate phase. If Wt > 1/3 
the scale-invariant amplitude increases and this occurrence can be appreciated both from 
the left and from the right plots of Fig. HI In the plot at the left of Fig. H] the amplitude is 
plotted as a function of the barotropic index of the intermediate phase. The lower limit for Wt 
coincides with —1/3 and it corresponds to a phase when the deceleration parameter vanishes. 
The upper limit on i^t is 1. In the case of constant barotropic index Cg^ = Wf In the limiting 
case Wt = I = and the speed of sound coincides with the speed of light. As argued in [23j, 
barotropic indices Wt > 1 would not be compatible with causality. Post-inflationary phases 
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Figure 4: The magnitude of the scale-invariant amphtude is illustrated as a function of Wt 
and as a function of (. 

stiffer than radiation (i.e. 1/3 < i^t < 1) can arise in different contexts (see, for instance, 
|24] ) where ( can be as small as 10~^^. If ti^t < 1/3 the conclusions are opposite: the scale 
invariant amplitude is always smaller than in the case when C = 1 when the transition to 
radiation is sudden. The results of Fig. IHshow that scale-invariant amplitudes ranging from 
0.01 nG to few nG cannot be theoretically excluded in the present framework. 

Let us now move to the case of blue spectral indices which are also the ones phenomeno- 
logically preferred at least heeding recent attempts aimed at constraining possible distortions 
of CMB anisotropies induced by large-scale magnetic fields. In this case the magnetic field 
intensity must be regularized and the commonly employed strategy is to use a Gaussian 
window function so that 

Bl = {B,{x,To)B,{x,ro)) ^ J ^(5.(^, ro)i?.(^, ro))e-'='^^ (4.6) 

where L = 2ti jki^. Equation (glD implies that Bl = ^B(2vr)i-"B r((nB-l)/2). In Fig. [5]the 
regularized amplitude is illustrated as a function of the spectral index and as a function of the 
barotropic index. The horizontal lines in both plots illustrate the requirements stemming 
from magnetogenesis demanding that the regularized amplitude must be, at least, larger 
than 

5L = 310=^e-^-*('^V^'nG, (4.7) 
Vpb/ 

where A^rot is the number of (effective) rotations performed by the galaxy since gravitational 
collapse; pa and pb denote, respectively, the matter density after and before gravitational 
collapse. The typical rotation period of a spiral galaxy is of the order of 3 x 10^ yrs which 
should be compared with 10^° yrs, i.e. the approximate age of the galaxy. The maximal 
number of rotations performed by the galaxy since its origin is of the order of A^rot ~ 30 (i.e. 
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the ratio of the two figures mentioned in the preceding sentence). This generous estimate 
refers to the situation where the kinetic energy of the plasma is transferred to the magnetic 
energy by means of an extremely efficient dynamo which amplifies the magnetic field by 
one efold for each rotation. The effective number of efolds is always smaller than 30 for 
various reasons. Typically it can happen that the dynamo quenches prematurely because 
some the higher wavenumbers of the magnetic field become critical (i.e. comparable with 
the kinetic energy of the plasma) before the smaller ones. Other sources of quenching have 
been recently discussed in the literature (see, for an introduction to this topic, section 4.2 of 
[25] and references therein). 

There is also another source of amplification of the primordial magnetic field and it 
has to do with compressional amplification. At the time of the gravitational collapse of 
the protogalaxy the conductivity of the plasma was sufficiently high to justify the neglect 
of nonlinear corrections in the equations expressing the conservation of the magnetic fiux 
and of the magnetic helicity. The conservation of the magnetic fiux implies that, during 
the gravitational collapse, the magnetic field should undergo compressional amplification, 
i.e. the same kind of mechanism which is believed to be the source of the large magnetic 
fields of the pulsars. Right before the gravitational collapse of the protogalaxy the mean 
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Figure 5: The regularized amplitude i?L is illustrated in the case of blue spectral indices. 

matter density can be taken to be of the order of the critical density, i.e. pb — Pcrit- After 
collapse the matter density pa of the galaxy is larger by, roughly, 6 or 7 orders of magnitude. 
Thus, after gravitational collapse, the magnetic field gets roughly amplified by a factor 
(Pb/Pa)^''^ — 10^ in spite of the possible occurrence of a strong or moderate dynamo action. 
It is worth mentioning that the scale-invariant amplitudes illustrated in Fig. H] satisfy the 
magnetogenesis requirements even in the absence of a strong dynamo action. 

In Fig. [5] the two horizontal lines in both plots illustrate the requirement expressed by 
Eq. f l4.7p in the case of two different numbers of effective rotations. In both plots of Fig. [5] 
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the regions which are above the horizontal hnes allow for a viable class of magnetogenesis 
scenarios. In the sudden reheating approximation (i.e. C = 1) the left plot of Fig. [5] shows 
that, for instance, ne < 1.35 to be compatible with the magnetogenesis requirement when 
A^rot = 20 (see thin oblique line in left plot of Fig. [5]). A relatively short delay in the 
beginning of the radiation epoch (i.e. ( ^ 10~^) allows for < 1.4 provided Wt > 1/3, 
i.e. provided the sound speed during the intermediate phase is larger than If, on 

the contrary, Wt < 1/3 the delay of the radiation results in a decrease of the magnetic 
field intensity and, ultimately, in a stronger bound on (i.e. n-Q < 1.3); this features are 
illustrated in the left plot of Fig. [5] with the dashed and with the dot-dashed oblique lines. 
In the right plot of Fig. Othe magnetogensis bounds are shown with the two horizontal lines 
(with naughts and diamonds). As written in the legend the delay of the radiation phase is 
assumed to be C ~ 10"'' and the magnetic field intensity is plotted as a function of Wt. 

Large-scale magnetic fields have been also observed in regular (or Abell) clusters and 
statistical determinations of intra-cluster magnetic fields are available since almost ten years 
[261 [271 [28] (see also [1] for a general perspective on cluster magnetism). These statistical 
determinations of cluster magnetic fields have been made possible thanks to the combination 
of Faraday rotation measurements with the determinations of the surface brightness of the 
clusters in the x-rays. The latter (satellite) determination allowed to infer the electron 
density and, ultimately, the strength of the magnetic field which turns out to be in the range 
of 500 nG. 

The typical scale of the gravitational collapse of a cluster is larger (roughly by one order 
of magnitude) than the scale of gravitational collapse of the protogalaxy. The mean mass 
density within the Abell radius ( ^ l.bh^^ Mpc) is roughly 10'^ larger than the critical 
density since clusters are formed from peaks in the density field. Moreover, clusters rotate 
much less than galaxies even if it is somehow hard to disentangle, observationally, the global 
(coherent) rotation of the cluster from the rotation curves of the constituent galaxies. By 
assuming, for instance, A^rot = 5, a density gradient of 10^ and 500 nG as final field, Eq. 
( 14. 7p demands and initial seed of the order 0.15 nG over a typical pivot scale /cl = ki,/ 10. 
This requirement, in the present context, holds when the spectral amplitude is quasi- fiat. 

In spite of the fact that it is always good to attempt accurate estimates, we must also 
admit that to have firmer predictions on the pre-decoupling magnetism it is necessary to 
model more carefully the evolutionary history of the conductivity between the end of inflation 
and the onset of the radiation-dominated epoch. This aspect is essential when assessing the 
viability of this type of scenario. 
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5 Concluding remarks 



In the last decade ground based as well as satellite observations probed with greater accuracy 
the nature of the magnetic field of the sun. The SOHO [29], TRACE [30], GONG [31] 
experiments |^ helped in deepening solar observations. In spite of these direct probes, it 
is fair to say that the 22-year magnetic cycle of the Sun cannot be claimed to be fully 
understood in terms of the so-called aio mechanism. The GONG observations give, for 
instance, the profile of the angular velocity of the sun not only on the surface but also in 
the interior of the sun. In spite of this valuable determination other effects may play a role 
so that it is fair to say that the situation is not clear [32j. 

The sun is the closest and better observed astrophysical object and still the features 
of dynamo amplification are under debate. In the case of large-scale magnetism we face a 
similar situation. It seems that we will probably never be able to know the initial conditions 
of the galactic dynamo as accurately as we pretend to understand the initial conditions of 
solar dynamos where, still, crucial puzzles remain. 

The hopes of clarifying this problem in the near future might not be so forlorn. Indeed, 
as repeatedly argued, CMB anisotropics and polarization [6l [25] are a powerful window on 
pre-decoupling physics and, consequently, on pre- decoupling magnetic fields. Waiting for 
more direct observational evidences it is plausible to speculate that a moderate dynamo ac- 
tion combined with compressional amplification could indeed bridge the regime of the initial 
conditions with the observed large-scale magnetic fields. Even in this simplified framework 
important theoretical problems persist and they have to do with the way the Universe be- 
comes a good conductor at the end of infiation. The early variation of the gauge couplings 
is a potential candidate for producing large-scale magnetic seeds for galaxies and clusters. 
However, more effort is certainly required especially in modeling the finite conductivity ef- 
fects and the self-consistent evolution of the whole system of equations. 



^ SOHO stands for Solar and Heliospheric Observatory; TRACE stands for TRansition Region and 
Coronal Explorer; GONG stands for Global Oscillation Network Group. 
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